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The metric perturbation induced by a particle in the Schwarzschild background is usually calculated in the Regge- 
Wheeler (RW) gauge, whereas the gravitational self-force is known to be given by the tail part of the metric perturba- 
tion in the harmonic gauge. Thus, to identify the gravitational self-force correctly in a specified gauge, it is necessary 
to find out a gauge transformation that connects these two gauges. This is called the gauge problem. As a direct 
approach to solve the gauge problem, we formulate a method to calculate the metric perturbation in the harmonic 
gauge on the Schwarzshild backgound. We apply the Fourier-harmonic expansion to the metric perturbation and 
reduce the problem to the gauge transformation of the Fourier-harmonic coefficients (radial functions) from the RW 
gauge to the harmonic gauge. We derive a set of decoupled radial equations for the gauge transformation. These 
equations are found to have a simple second-order form for the odd parity part and the forms of spin s = and 1 
Teukolsky equations for the even parity part. As a by-product, we correct typos in Zerilli's paper and present a set 
of corrected equations in Appendix M. 



A compact object of solar-mass size orbiting a supermassive black hole is one of promising candidates for the source 
of gravitational waves. Since the internal structure of such a compact object may be neglected in this situation, we 
may adopt the black hole perturbation approach with the compact object being regarded as a point particle. In the 
black hole perturbation approach, we consider the metric perturbation induced by a point particle of mass fi orbiting 
a black hole of mass M, where /i ^ M. At the lowest-order in the mass ratio (fi/M), the motion of the particle 
follows a geodesic of the background spacetime. In the next order, however, the particle moves no longer along a 
geodesic of the background because of its interaction with the self-field. 

Although this deviation from a background geodesic is small for n/M ^ 1 at each instant of time, after a large 
lapse of time, it accumulates to become non-negligible. For example, a circular orbit will not remain circular but 
becomes a spiral-in orbit and the orbit eventually plunges into the black hole. 

If the time scale of the orbital evolution due to the self-force is sufficiently long compared to the characteristic orbital 
time, we may adopt the so-called adiabatic approximation in which the orbit is assumed to be instantaneously geodesic 
with the constants of motion changing very slowly with time. In the Schwarzschild background case, we may assume 
the orbit to lie on the equatorial plane and the geodesic motion is determined by the energy and (the z-component 
of) angular momentum of the particle. In this case, the time variation of the energy and angular momentum can be 
determined from the energy and angular momentum emitted to infinity and absorbed into the black hole horizon by 
using the conservation law. 

However, there are cases when the adiabatic approximation breaks down. For example, in the case of an extremely 
eccentric orbit or an orbit close to the inner- most stable circular orbit, the orbital evolution will not be adiabatic 
because the stability of the orbit is strongly affected by an infinitesimally small reaction force. Furthermore, in the 
Kerr background, there is an additional constant of motion, known as the Carter constant. Intuitively, it describes 
the total orbital angular momentum, but unlike the case of spherical symmetry, it has nothing to do with the Killing 
vector field of the Kerr geometry. The lack of its relation to the Killing vector makes us impossible to evaluate the 
time change of the Carter constant from the gravitational waves emitted to infinity and to event horizon, even in the 
case when the adiabatic approximation is valid. Thus it is in any case necessary to derive the self-force of a particle 
explicitly. 

The gravitational self-force F^^ is formally given as 
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where {z"(r)} represents the orbit with r being the proper time measured in the background geometry and F"^ 
the connection of the background. The self-force arizes from the metric perturbation h^^^ induced by the particle: 



and it is expressed as 



where Pa^ = 6a^ + UaU^ , hap — hap — \gaph and u°' — dz'^ /dr. 

The metric perturbation diverges at the location of the particle and so does the self-force. Thus the above formal 
expression is in fact meaningless. Fortunately, however, it is known that the metric perturbation in the vicinity of the 
orbit can be divided into two parts under the harmonic gauge condition; the direct part which has support only on 
the past light-cone emanating from the field point x^^ and the tail part which has support inside the past light-cone, 
and the physical self-force is given by the tail part of the metric perturbation which is regular as we let the field point 
coincide with a point on the orbit; x^^ — > 2;^(r) j^, ^. It must be noted that the direct part can be evaluated by 
local analysis, i.e., only with the knowledge of local geometrical quantities. Therefore the physical self- force can be 
calculated as 

lim Fa[h'^^{x)]^ lim {Fa[h{x)] ~ Fa[h'^" (x)]) . 

Furthermore, it has been revealed recently by Detweiler and Whiting |^ that the above devision of the metric can be 
slightly modified so that the new direct part, called the S part, satisfies the same Einstein equations as the full metric 
perturbation does, and the new tail part, called the R part, satisfies the source-free Einstein equations, and that the 
R part gives the identical, regular self-force as the tail part does. The important point is that the S part can be still 
evaluated locally near the orbit without knowing the global solution. 

When we perform this subtraction, we must evaluate the full self-force and the direct part under the same gauge 
condition. But the direct part is, by definition, defined only in the harmonic gauge. On the other hand, the full metric 
perturbation is directly obtainable only by the Regge-Wheeler-Zerilli or Teukolsky formalism |^, |^. Therefore 
one must find a gauge transformation that brings both the full metric perturbation and the direct part of the metric 
perturbation to those in the same gauge. This is called the gauge problem Q. 

For the direct part, methods to obtain it under the harmonic gauge condition were proposed & |[ |ll|. However, 
it seems extremely difficult to solve the metric perturbation under the harmonic gauge because the metric components 
couple to each other in a complicated way. This is one of the reasons why the gauge problem is difficult to solve. 

Recently, Barack and Ori [13] gave a useful insight into the gauge problem. They proposed an intermediate 
gauge approach in which only the direct part of the metric in the harmonic gauge is subtracted from the full metric 
perturbation in the RW gauge. They then argued that the gauge-dependence of the self-force is unimportant when 
averaged over a sufficiently long lapse of time. Using this approach, the gravitational self-force for an orbit plunging 
into a Schwarzschild black hole was calculated by Barack and Lousto |13|. But they also pointed out that the RW 
gauge is singular in the sense that the resulting self-force will still have adirection-dependent limit for general orbits. 
The situation becomes worse in the Kerr background where the only known gauge in which the metric perturbation 
can be evaluated is the radiation gauge but the metric perturbation becomes ill-defined in the neighborhood of 
the particle, i.e., the Einstein equations are not satisfied there [T^ . 

In this paper, as a direct approach to the gauge problem, we consider a formalism to calculate the metric perturbation 
in the harmonic gauge. We focus on the Schwarzschild background. Instead of directly solving the metric perturbation 
in the harmonic gauge, we consider the metric perturbation in the RW gauge first, and then transform it to the one in 
the harmonic gauge. Namely, we derive a set of equations for gauge functions that transform the metric perturbation 
in the RW gauge to the one in the harmonic gauge. 

The paper is organized as follows. In Sec. II, we formulate the gauge transformation from the RW gauge to the 
harmonic gauge. First we decompose the gauge transformation generators into the Fourier-harmonics components. 
Then the generators are devided into three parts; the odd parity part and the even parity part which is further devided 
into scalar and di verg ence free parts. By the above procedure, we find a set of decoupled equations for the gauge 
functions. In Sec. IIF we summarize our formulation and discuss remaining issues. In Appendix we recapitulate 



the equations for the Regge-Wheeler-Zerilli formalism by correcting typos in Zerilli's paper 

II. FORMULATION 

We consider a metric perturbation h^i, in the Schwarzschild background, 

gy.^dx^dx'' ^-!{r)dt' + ]{r)-^dr'^ ^r^{dd^ ^^\^^ed<l?); /(r) 1 - ^ . (2.1) 
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We express the gauge transformation from the RW gauge to the harmonic gauge as 

^RW ^ ~ ^RW ^ ^'^ ' (2-2) 

lRW ^ i,H _ , RW _<:_/: (2 -i) 

where the suffix RW sta nds for the RW gauge and H for the harmonic gauge. 

Substituting Eq. (2.3) into the harmonic gauge condition h^^y''^ = 0, we obtain the equations for f : 

f \i' _ IRWi/ _ i,RWiy _^ 10 A\ 

Using the static and spherical symmetry of the background, we perform the Fourier-harmonic expansion of the above 
equation and consider the equations for the expansion coefficients for and . We use the tensor harmonics 
introduced by ZeriUi Q which are recapitulated in Appendix Then according to the property under the parity 
transformation (0, (/i) ^ (vr — 0, + tt), we devide into the odd and even parity parts, 

^ ^(odd) + ^fcvcn) ' (2.5) 

and the even part is further decomposed into the scalar and divergence-free part, 

?reve„)-^^"+^rv)' (2.6) 

where ■, =0. 



A. Odd part 



First, we consider the odd part which has the odd parity (—1)^"''^ under the parity transformation. The gauge 
transformation generators and the metric perturbation are given in the Fourier-harmonic expanded form as 



f(odd) 



, (odd) 
iifiv 



j dio e 

/J™ ^ 



' smd 



(2.7) 



hQlrnuj{r)Cf.^ + ih 



V(^-l)(^ + 2) 



2r 



(2.8) 



where c^^^j^,C£m^i^, dim^u are tensor harmonics with odd parity 1^. By substituting Eqs. (2.7) and (2.8) into Eq. ( ^.4| ), 
we obtain the equation for Aimuiir)- 



16 i TTr 



; variable, 

stress-energy tensor, and the differential operator C is defined as 



(2.9) 



where R^m!^^^''') ^^'^ Regge- Wheeler gauge invariant variable, DimLj{r) is the Fourier- harmonic coefficient of the 





ar 


( 

H/(r) 


1 


( 




/(r) dr* ^ 


^ \f{r) 


r2 



(2.10) 



where r* = r + 2Mlog(r/2M — 1) and ' — d/dr. The form of the differential operator C is slightly different from the 
radial part of the scalar d'Alembertian, but Eq. (2.9) may be solved by the standard Green function method. The 
homogeneous solutions to construct the Green function can be obtained by applying the method developed by Mano 
et al. . Here we note that the retarded causal boundary condition should be imposed for the Green function. 
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B. Even scalar part 



Next we consider the even scalar part of which has the even parity (—1)^ and is expressed by a gradient of a 
scalar function ^. It is expressed as 



(2.11) 



The gauge equation for ^j^^ is derived from ^jf-l'" = J^^\ which, because of the vanishing of the background Ricci 
tensor, gives 



where the source term J^^^ is determined from the equation, 



(2.12) 



(2.13) 



The Fourier- harmonic expanded form of the above equation is derived as follows. First, we take the divergence of the 
metric perturbation under the RW gauge; 



+ 



V^(£ + l)(^-l)(£ + 2)/2 



^imu;(r)(e3)^ > , 



(2.14) 



where (ct)^ = (—1,0,0,0), (e^)^ = (0,1,0,0) and (63)^ = (0, 0, SeY^m, 9,^Yem). Then we take the divergence of the 
above once more to obtain 



-£(^)F^„„(r)+87r 



(2.15) 



where Hemw{i') is the trace of the metric perturbation given by 



J^emu{r) — J^emuV) ■ 



(2.16) 



The d'Alembertian of J^^^ is expanded as 



e— *y,„((^,0)/:(^)j(l(r) 



(2.17) 



where ji^^^ is the Fourier-harmonic coefficient of J^^^ and C'^^^ is the radial part of the scalar d'Alembertian, 



,(.).,(,)^ + ;V)^^( 



u;2 £(£+1) fir) 



_J. (f_ / g;^ £{£+!) 2M 

f{r)dr*^ \f{r) 



(2.18) 
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Hence Eq. ( 2.13| ) becomes 



■muj 



Once the source term is obtained, we obtain from Eq. ( 2.12 ) the equation for £,£rnuj{r) as 



(2.19) 



(2.20) 



Note that this equation as well as Eq. (2.19) are just the radial part of the scalar d'Alembertian, or equivalently the 
spin s = Teukolsky equation, hence can be also solved by the Green function method. The detail analysis of the 
homogeneous solutions which satisfy the s — Teukolsky equation is discussed in H. 



C. Even vector part 



Since the even vector part ^^^^ satisfies the divergence free condition, ^^^^'^ = 0, this part has two degrees of 
freedom. Therefore ^^^^ is expressed in terms of two independent radial functions. 



-Motmu,{r)Ye„^{0, (j)){et)^ + -Afiw(r)Yc,„(6', (?!))(e^)^ 



+ . <! ^r^Moimu^ir) + dr {rf{r)Mum^{r)) } (eg)^ . (2.21) 

To solve the gauge transformation equation ( ^.4| ) for this part, we introduce an auxiliary field F^i, := £}j^^ — ^jx^l- 
Then we have 



p _ _ t(v) . t(v) ^ 7RW:i/ _ j(s) 

^ fJ-t^ ^ ^/i I ^/i — "-III/ "^lU I 



(2.22) 



where we have again used the fact that the background Ricci tensor vanishes. Note that the Fourier-harmonic 
expansion of h^'^''^ is given by Eq. (^.14 ). The above equation is the Maxwell equation, so it can be solved by using 
the Teukolsky formalism for the electromagnetic field (s = ±1). Namely, we introduce the following variables: 



where * denotes the complex conjugate, and , n**, m** are the Kinnersley null tetrad defined by 



1 



fir) 



1,0,0 , -,- 



1 f{r) 



0,0 



1 



^/2r 



0,0,1,- 



sin 9 



The variables (pQ and 02 satisfy the s = ±1 Teukolsky equation, 

^(Tcuk) ^ _47rr2T, , 

where = 4>o and ~ (1)2 / , and is the Teukolsky differential operator defined as 



(2.23) 



(2.24) 



(2.25) 



£(Touk) 



-dl + 2s 



M 
W) 



r]dt + {r^f{r)) ' dr {r^ f (r))'''' 8, 



: ^ 1 f ■ aa \ : ^ a2 , 2iscos6' 
sm sm 6 ^ sm 9 



The source terms Tj , which are calculated from jj^'^ , are 



Ti 



T_i = - 



/2r 
1 



V2r 



do- 



sin 9 



f{r] 



— s(s cot^ ^ ^ 1) 















^ fir) 



(2.26) 

(2.27) 
(2.28) 
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The solution of this Teukolsky equation is also analyzed in . 

Once we obtain (/jq and (1)2, the two radial functions for the gauge transformation are calculated from Eq. (2.23) 
which reads 









£m 


—iuj < 


—iur 


I f{r) 


Jd. 






tin 


+iu; ■ 




I f{r) 



lYlnAOA) 



MoernUr) + ■^{rf{r)Mu„r^{r)) \ + f{r) 



d 



-ZUJT d 



iojt -lYi,n{0,( 



MoimUr) + ^(r/(r-)Mw™c.(r)) } + f{r) — 



HOT d 



(2.29) 



(2.30) 



where sYimi^T't') sse the spin- weighted spherical harmonics. Performing the Fourier and spin-weighted spherical 
harmonic expansion for 0o and 02, 



ho = dLu'^4>oi„rujir)e "^\Yt„, 
h ^ J duY,4>2ir,rUr)e"^'~iYi 



we obtain the equations, 



2 ~ 



rf{r) 

v/2/^(^+l) 
rf{r) 



l{e+l) uj^r\ , , d , ^ ^ 

-JTT MQim^[r) - tLU — [rf(r)Mumuj[r)) 

r Jv)/ dr 

i{e+l)fir 



Mu„iuj{r) 



'fir)-f^irfir)MumUr))+^^f{r)^ . , 
dr'^ dr \j{r) 



Moimojir) 



(2.31) 
(2.32) 

(2.33) 
(2.34) 



We can eliminate Murnu>{r) from Eqs. (p.33|) and (2.34) to obtain a decoupled equation for Ma^rnu>{r) 



fir)'-^^oimuir)+f{r) 



2 + 2/(r) , . \ d ~ 

+ fir) + luj —cpoemu^ir) 

r I dr 



, rf'ir) -{l-m + 2) fir) + ^ur{l + 2/(r)) - x^'l 
H 5 (PQemLu[r) + 2i{r)-—(j)2tmuj{r) 



2 + 2/(r) ■ W I I . 

r J dr 



dr'^ 

rf (r) + (£ - 1) (£ + 2) + Situr ~ 



(t)2irnLo{r) 



,(2.35) 



47r 2r2f_i,„„(r) + /(r)f«„^(r) 



2rV(r) 

f( \c , fUW*^! /\ {'^ + -i^^r){rf'{r)+iujr)~ 
dr 



+ -iiruj - 1 - 7fir))-^4>2imu.ir) + 2 
r dr 



ui"^ + iui/r 5/(r) — 1 



fir) 



+ 



02£ma, (r) 



(2.36) 



wher e T^f rr,.,,)(r) is the Fourier-harmonic coefficient of the source term in the Teukolsky equation ( p.25| ). Interestingly, 
Eq. ( 2.36 ) has the same form as the s = Teukolsky equation. Once we obtain Moimujir) by solving the above 
equation, Mumuiir) is derived from 



Miimu>ir) 



ujy/2£{e+l) 



d ~ , . 1 + iujr ~ r. d 2 / \ n /(^) ~ 7 I \ 

}\r)-r9WTaui\r) ^ q>oimu\r) + 2 — (p2imu\r) -I- 2 — - — (j)2e7nuj{r) 

dr r dr rjyr) 



( d 



—MiMrnu,{r) - 



Moini^ (r) 



V dr 



(2.37) 



which also follows from Eqs. (2.33) and (2.34). 
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III. SUMMARY AND DISCUSSION 



In this paper, to solve the gauge problem of the gravitational self-force, we have considered the gauge transformation 
from the Regge- Wheeler gauge to the harmonic gauge and have presented a formalism to obtain the infinitesimal 
displacement vector of this transformation, First, we have performed the Fourier- harmonic expansion of and 
divided it into the odd and even parity parts. The odd part has only one degree of freedom and it turns out that 
the gauge transformation can be found by solving a single second-order differential equation for the radial function. 
As for the even parity part, we have further divided it into scalar and vector parts where the scalar part is given 
by the gradient of a scalar function and the vector part is divergence-free. The scalar part has by definition only 
one degree of freedom, and we have found that it can be obtained by solving two second-order differential equations 
consecutively. These two equations are found to be identical to the s — Teukolsky equation. The vector part has 
two degrees of freedom, and the gauge transformation equations give equations that are coupled in a complicated way. 
However, by introducing two auxiliary variables which satisfy the s = ±1 Teukolsky equations, we have succeeded in 
deriving a decoupled second-order equation for one of the gauge functions with the source term given by the auxiliary 
variables. Interestingly, this second-order equation has the same form as the s = Teukolsky equation. The other 
gauge function is then simply given by applying a differential operator to the first. 

Since all the equations to be solved have the form analogous to or equal to the Regge- Wheeler equation, we can 
derive analytic expressions for their homogeneous solutions by using the Mano-Suzuki-Takasugi method [|l4| , and 
construct the Green function from these homogeneous solutions. So we conclude that the gauge transformation can 
be solved by using the Green function method, and we can construct the metric perturbation in the harmonic gauge. 
In practice, however, it may not be easy to solve for the gauge transformation since it involves products of Green 
functions with double integrals. Derivation of the gauge transformation functions in a closed, practically tractable 
form is left for future study. 

Another approach to the gauge problem is to consider the self-force in a gauge different from the harmonic gauge, 
similar to (but very different in principle from) the intermediate gauge approach proposed by Barack and Ori [ p2[ . 
Here the recent result by Detweiler and Whiting j|] becomes crucial. Their observation that the S part and theR 
part play the identical roles as the direct part and the tail part, respectively, and that the S part satisfies the same 
inhomogeneous Einstein equations as the full metric perturbation enables us to define the S part and the R part 
of the metric perturbation unambiguously in an arbitrary gauge as long as the gauge condition is consistent with 
the Einstein equations. For example, given the S part of the metric perturbation in the harmonic gauge, one can 
perform the gauge transformation of it to the RW gauge and the resulting metric perturbation which satisfies the 
Einstein equations can be identified as the S part of the metric perturbation in the RW gauge. Then, after solving 
the Regge- Wheeler- Zerilli equations to obtain the full metric perturbation, it is straightforward to derive the R part 
of the metric perturbation in the RW gauge The calculation of the self-force in the RW gauge in this manner is 
in progress |16|| . 

Finally, we comment on the self-force in the case of the Kerr background. In the Schwarzschild case, it was possible 
to use the Regge- Wheeler- Zerilli formalism to obtain the metric perturbation in the RW gauge. However, in the Kerr 
case, there is no known gauge in which the full metric perturbation can be calculated. The Chrzanowski method 
Q based on the Teukolsky formalism can give the metric perturbation in the (ingoing or outgoing) radiation gauge, 
but only outside the range of radial coordinates the orbit resides in. One possible way to circumvent this difficulty 
is to consider first the regularization of the Weyl scalar '^4. Given an orbit, ^"4 can be calculated by the Teukolsky 
formalism, and the S part of it, can be calculated from the S part of the metric perturbation in the harmonic 
gauge, hl'^, 

*! - ^4hl^^] , (3.1) 

where is the operator to derive the Wely scalar from a given metric perturbation. Then the R part of ^'4 can be 
derived by subtracting the S part from the Weyl scalar, 

= ^-4 - ^-f . (3.2) 

Now ^I'J satisfies the homogeneous Teukolsky equation. Hence using the Chrzanowski method, we may construct the 
R part of the metric perturbation in the radiation gauge and derive the self-force. Since this procedure involves many 
derivative operations, the metric perturbation h^'^ has to be evaluated to with a sufficiently high accuracy which 
may be practically a difficult task, if not impossible. Feasibility of this method should surely be investigated. 
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APPENDIX A: THE FIELD EQUATION ON THE REGGE- WHEELER GAUGE 

In this appendix, we recapitulate the equations for the Regge- Wheeler- ZeriUi formalism. In doing so, we correct 
some minor errors in Zerilli's paper Here an equation number given as (Z:l) denotes the equation (1) in Zerilli's 
paper for comparison, and a label [CRTD] to an equation means it is corrected. 

We consider the linearized Einstein equations for the perturbed metric, 

where g^^ is the background metric. Then the Einstein tensor and the stress-energy tensor up to the linear order can 
be expanded as 



Gfj,^[gfj.u] = Gfj^uig^iu] + SGfj,u[hf_cu] + 0{h^), 



(Al) 
(A2) 



where = hf^a'" and 



(A3) 



When the background is Ricci flat, R^^} — 0, the above equation is rewritten as 



(A4) 



We apply the above to the case of the Schwarzschild background, and expand ft-^i/ ((Z:D2a) and (Z:D2b)) and ST^i, 
in tensor harmonics. 



i 
r 



ST 



v/2^(£+l)(£ -!)(£ + 2) ^ 
2r 



[CRTD] 



[4(0)^(0) .(1) (1) . 4_ r(o)j,(o) I R h ^n(°'>J°'>^n n 



+ D£mdim + Gn Qirn + Fhnfl 



£m 



(A5) 



(A6) 



where we use /ig^^^ and h'f^^^ for the even part coefficients instead of /iq™^ and /Iq^, respectively, in Zerilli's paper, 

and the coefficient instead of the Zerilli's notation Gim for the energy-momentum tensor, and a^^^^, agrn , • ■ ■ are 
the ten tensor harmonics (Z:A2a-j) defined as 



(0) 



(1) - 






000/ 



i^/V2) 



( Ytr,^ 

Syva 



V 



(A7) 



(A8) 
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TABLE I: Stress-Energy Tensor in terms of Tensor Harmonics 



Description 



Dependence of "driving term" on r and t 

-2 



Tensor harmonic 



Even 



Even 



(„) / 2M\ 



^5{r - R{t))YUm) 



Even 



A 



(1) 



^y^if.-y-^r-^Sir ~ R{t))Yl^{n{t)) [CRTD] 



Even 



B 



(0) 



U{£ + l)]~'/'^/i7 (l - ^) r-'S{r - R(t))dY;U^{t)) / dt [CRTD] 



Even 



gftn = [\m + l)]~'^V7(r - 2M)-^^5 {r - R{t))dYtUn{t))/dt 



bim{9,(t)) 



Odd 



[^£{e + 1)1-^/^7 (l - ^) r-^<5(r - 7?(t)) 

1 ay,:„de ^i^Q.^y/™ di>" 



sin 6 9$ dt 



Odd 



1 rJR 

Q„n = + l)]-'''i^^^^{r - 2M)-'S{r - R{t)) 



1 dYl^dQ . ^dYl^d^ 
sm e 



sin 6 9$ dt 



d& dt 



[CRTD] 



Odd 



Dfm ~ 



' + l){e.-l){e. + 2)]-^/'i/i7<5(r - R{t)) 



'dO 2 ■ 2^fd$ 2 



sm B dt dt ; 



Even 



/d$ de 1 



+ 1)(£ - + 2)]-'/V7'5(r - R{t)) 

wu^iit)] 



Even 



G'';i = ^5{r~R{t)) 



,de 2 , . 2^,d$ 2 



Substituting Eqs. (A5) and ( |A6| ) into Eq. (A4), we obtain the field equations for each harmonic mode. For the odd 
part which has the odd parity (—1)^+^, in the RW gauge in which /12 — 0, the following three equations (Z:C6a-c) are 
derived. 



d^ho 




2dhi 


Qj.2 


dtdr 


r~dt 




d'ho _^ 


2dho 


dt^ 


dtdr 


r dt 






d 








dr 


[(■ 



4A/ e{£ + 1) 

2 







r-2M 



-hi 



r dho 
r-2M~dt 



^l{£+l)/2 r-2M ^'^ t^^^^]' 
87ri(r- 2Af) 



^^(^+l)(^-l)(^ + 2)/2 



(A22) 

CRTD] , (A23) 
D,™ [CRTD]. (A24) 



Fo the even part which has the even parity (—1)^, we have seven equations (Z;C7a-g) in the RW gauge in which 
'-0 



hi'^ = h['^ = G = 0. 



r I dr'' r \ '"/V r J dr 



1 / 2M\ dH2 1 A 2M\ 

dr 



7' / 



(i^2 - K) 



r I 



d_ 
dt 



{H2 + K) 
M 



^ + \k-H2)-- — - 

dr r r[r — 2M) 



2r2 



87r42 [CRTD], 



K 



Hi = 



-4V27ri42 [CRTD] , 



(A25) 



(A26) 
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^ d^K r-M dK 2 dH^ 
r-2M) dt^ ~ r{r - 2M) dr ~ r-2M "df 
I ^dHo 1 



2r(r - 2M) 



(i^-ifo) = -SnAem [CRTD], (A27) 



9r 



2M 



Hi 



2M. 

^ , Hr, - AC 1 -I- 

9t \ r J dr 



r 

r d^K / 2M\ S^/iT 2 / M\ SiC 



r — 2M dt^ dtdr \ r J dr^ 
^2{r-M)dHi 1/ M\dH2 r + M dHo 



r{r — 2M) dt r \ r J dr r"^ dr 



+ ^1^[H,-H,) = 8^/27^G(^ [CRTD] 



(A30) 



H0-H2 ^ 8nr^F,^ [CRTD]. (A31) 

2 ^/e{e +!){£ -!){£ + 2) /2 ^ i y ' 

We now consider the Fourier transform of the above field equations. The Fourier coefficients are defined, for 
example, as 

iuJt 



hoimUr)= dthoem{t,r)e''^- (A32) 
J —00 

Then we derive the Regge-Wheeler-Zerilli equations and construct the metric perturbation under the RW gauge 
condition. 

For the odd part, a new radial function R''^^\r) is introduced, in terms of which the two radial functions /lo^mw 
and humoj for the metric perturbation are expressed as 

hi,^. = (A33) 
id n(odd)x 87rr(r — 2M) 



The new radial function satisfies the Regge- Wheeler equation (Z:ll), 

+ [.^-vr\r)]Rrr 

SiTi r-2M 



[iZ(€+l)(£-l)(€ + 2)]i/2 r2 

X - ^)Demu.] + {r- 2M)[(^ - !){£ + 2)]V2q,^^^ [CRTD] , (A35) 

where r* = r + 2Mlog(r/2M - 1) and 

>'r'(-) = (i-^)(^^-^). (A3e) 
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For the even part, a new radial function R^^^^ (r) is introduced, in terms of which the four radial functions (Z: 13-16) 
are expressed as 



A(A + l)r^ + 3AMr + 6Af ^ (even) r - 2M dR^ZZ^ 
r2(Ar + 3M) ^"'^ r dr 



r{r - 2M) 



Clin 



i{r - 2Mf 



H 



\r + 3M """" r(Ar + 3M) 
Ar2 - 3AMr - SM"^ 



{r ~ 2M){Xr + U'l) 



R 



C2emu [CRTD] , 



imuj 



Emu 



dr 



Xr + 3M 



Lor{r - 2M) ~ 



Xr + 3M 



C2tmu. [CRTD] 



Xr{r - 2M) - uj^r'^ + M{r - 3M) M(A + 1) - uj\^ 



(r ~ 2M){Xr + 3M) 
where we have introduced the symbol A for 



16nr^[^i{£ +!){£- 1){£ + 2)]-^''^Fi. 



iujr{Xr + 3M) 
[CRTD] 



B 



Imu) 



1 



X = -{t-l){l + 2), 
and the local source terms ((Z:17), (Z:20) and (Z:21)) by 



(A37) 
(A38) 

[CRTD] , (A39) 
(A40) 

(A41) 



Climui 
C2lmui 



Xr + 3M 
8^ ^^(1) 



1 ^ 
— J 

r 



V2(jJ 

iuj r - 2M 
Wirir^ A 



B. 



£muj 



(0) _ 



r-2M 



-Be 



2M'2 



-£{£ + - 1)(£ + 2)]-i/2j^,„^ [CRTD] , 



(A42) 
(A43) 

(A44) 



We note that the above radial functions for the metric perturbation have the local source terms which have the 
i5-function behavior at the particle location. The new radial function obeys the wave equation, 



d2 D (even) 



2 -rr(oVOn)/ N 

- (r) 



R 



(even) 



£mLU 5 



where 



Tr(even)^ n 



2M\ 2A2(A + l)r3 + GX^Mr^ + ISXAPr + 18M^ 
~J r3(Ar + 3M)2 



(A45) 



(A46) 



and the source term is 



+i 



-~2M d 
r dr 
(r-2M)2 



(r - 2Mf 



-.c 



limuj 



r{Xr + 3M)2 



r{Xr + 3M) \r -2M 

A(A + l)r2 + 3AMr + GM^ ~ 
^ C' 



C: 



2tn 



2ln 



+ i- 



. Xr^ - 3XMr - 3M^ ~ 



2M 



limuj 



(A47) 



The above equation ( A45) is called the Zerilli equation. The Zerilli equation can be transformed to the Regge- Wheeler 
equation by the Chandrasekhar transformation |l^ . So we may focus only on the Regge- Wheeler equation if desired. 
The Regge- Wheeler homogeneous solutions are discussed in detail by Mano et al. |Q. Using their method, one 
can construct the retarded Green function to solve t he in homo geneo us Regge- Wheeler equation. Th en the me tric 
perturbation in the RW gauge is obtained from Eqs. ( A33| ) and (A34) for the odd part and from Eqs. (A37- A40| ) for 
the even part. 



[1] Y. Mino, M. Sasaki and T. Tanaka, Phys. Rev. D 55, 3457 (1997) [ MXiv:gr-qc/960601^ ] 



13 



[18] 



T. C. Quinn and R. M. Wald, P hvs. Rev. D 60. 0640 09 (1999) | |arXiv:gr-qc/9903014 |. 

S. Detweiler and B. F. Whiting, |arXiv:gr-qc/0202086| . 

P. L. Chrzanowski, Phys. Rev. D 11, 2042 (1975). 

T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). 

F. J. Zerilh, Phys. Rev. D 2, 2141 (1970). 

S. A. Teukolsky, Astrophys. J. 185, 635 (1973). 

L. Barack and A. Ori Phys. Rev. D 61, 061502 (2000) | |arXiv:gr-qc/9912010| , L. Barack, Phys. Rev. D 62, 084027 (2000) 
|irXiv:gr-qc/0005042[ . 

Y. Mino, H. Nakano and M. Sasaki, |arXiv:gr-qc/0111074 

L. Barack, Y. Mino, H Nakano, A. Ori and M. Sasaki, Phys. Rev. Lett. 88, 091101 (2002) | |arXiv:gr-qc/0111001 |. 
Barack and A. Ori, |arXiv:gr-qc/0204094 

Barack and A. Ori, Phys. Rev. D 64, 124003 (2 001) | |arXiv:gr-qc/0107056t . 
Barack and C. O. Lousto, arXiv:gr-qc/0205043 



Mano, H. Suzuki and E. Takasugi. Prog. The or. Phvs. 95. 1079 (1 996') | }a.rXiv:gr-qc/960302C| ], S. Mano, H. Suzuki and 
E. Takasugi, Prog. Theor. Phys. 96, 549 (1996) larXiv:gr-qc/9605057||. 

Y. Mino, talk given at the 5th Capra Ranch meeting, |http: / /cgwp. gravity. psu.edu /events/Capra5 /slides / mino.ppt 
H. Nakano, N. Sago and M. Sasaki, in preparation. 

S. Chandrasekhar, Proc. R. Soc. (London) A343, 289 (1975); Mathematical Theory of Black Holes (Oxford University 
Press, 1983), 

M. Sasaki and T. Nakamura, Phys. Lett. A 89, 68 (1982). 

Furthermore, the gravitational self-force is, because of the equivalence principle, a gauge- variant notion. To give a genuinely 
physical meaning to it, one must solve the second order metric perturbation completely. This is however beyond the scope 
of the present paper. 



